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A Dynamical  Model  for  the  Onset  of 
Magnetospheric  Substorms 

I 

1.  INTRODUCTION 

■ 

[ 

The  purpose  of  this  paper  is  to  develop  a dynamical  model  for  the  onset  of 
magnetospheric  substorms  based  on  experimental  observations  by  Pytte  et  al.  * 

We  require  this  model  to  contain:  (1)  a "triggering"  mechanism  that  is  related 
to  the  compression  of  the  plasma  sheet;  (2)  a merging  process  which  is  governed 
by  magnetic  (current)  fluctuations;  and  (3)  a mechanism  that  explains  multiple 
onsets. 

Present  steady-state  merging  theories  fall  into  two  general  categories: 

(1)  self-consistency  conditions  between  the  convected  particles  and  the  crosstail 

2 3 • 

current  (Alfven,  and  Dessler  ),  and  (2)  dissipation  models  for  the  conversion  of 

4 5 6 

magnetic  energy  into  joule  heating  (Sweet,  Parker,  and  Coppi  and  Friedland  ). 

An  excellent  review  article  by  Vasyliunas  also  treats  the  standing  shock  theory  of 

g 

Petschek. 

Approach:  The  behavior  of  the  plasma  sheet  is  governed  by  the  relative  values 
of  the  mean  collision  time,  the  proton  gyro-period  and  the  characteristic  transit 
time  of  the  J X B force.  The  latter  will  be  explained  in  detail  later.  We  first 
assume  and  then  show,  that  if  the  mean  proton  collision  time  is  sufficiently  long 


(Received  for  publication  11  December  1978) 

Due  to  the  number  of  references  to  be  included  as  footnotes  on  this  page,  the 
reader  is  referred  to  the  list  of  references,  page  31. 
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that  a finite  number  of  protons  can  travel  across  the  plasma  sheet  without  being 

diverted  by  collisions.  These  protons  are  simply  responding  to  the  J X B force. 

g 

With  the  magnetic  field  measurements  of  Bowling  and  Wolf,  the  J X B equa- 
tion of  force  turns  out  to  be  that  of  a harmonic  oscillator.  It  is  shown  that  in 

2 15 

order  for  the  harmonic  component  to  be  significant  that  n i < 8.  1 X 10  protons/ 

3 

cm.  For  n ~ 0.  3 protons /cm  , this  condition  is  satisfied  if  the  north-south 
thickness  of  the  plasma  sheet,  i,  is  less  than  1634  km.  Thus  compression  of  the 
plasma  sheet  by  the  tail  lobes  enhances  the  harmonic  component.  This  in  turn  is 
identified  as  a time-dependent  "triggering"  mechanism.  An  individual  particle 
treatment  given  in  Appendix  A leads  to  a similar  conclusion. 

In  Section  3,  the  effects  of  a localized  interruption  of  the  cross-tail  current 
are  analyzed.  If  the  interruption  occurs  near  the  neutral  sheet,  then  symmetric 
magnetic  field  line  deformations  will  occur  at  the  plasma  sheet -tail  lobe  interfaces. 
If  the  harmonic  oscillator  triggering  mechanism  is  present,  the  deformed  magnetic 
field  lines  may  merge.  It  is  found  that  one  of  the  necessary  conditions  for  the 
process  is  related  to  and  thus  controlled  by  the  daytime  merging  rate. 

In  Section  4,  we  develop  a simple  dynamical  model  for  multiple  onsets.  Force 
constants  for  the  tail  lobes  are  derived  by  treating  transverse  compressional  waves 
as  being  equivalent  to  a density  of  harmonic  oscillators.  These  force  constants, 
coupled  to  that  of  the  simple  harmonic  component  in  the  plasma  sheet  (Section  2), 
represent  a three-oscillator  system.  The  eigenmodes  of  this  system  are  derived 
in  the  standard  fashion.  One  of  the  four  possible  eigenmodes  corresponds  to  a 
displacement  of  both  tail  lobes  towards  the  neutral  sheet.  It  is  this  mode  one 
expects  to  be  excited  during  magnetic  substorms.  The  corresponding  eigenperiod 
is  a sensitive  measure  of  the  simple  harmonic  oscillator  component  in  the  plasma 
sheet  (Section  2)  and  of  the  mean  proton  collision  frequency. 

The  physical  picture  of  a substorm  according  to  this  model  is  as  follows: 
Simultaneous  displacement  of  both  tail  lobes  towards  the  neutral  sheet  compresses 
the  plasma  sheet.  When  this  compression  is  sufficiently  intense,  a harmonic  com- 
ponent can  then  be  maintained  across  the  neutral  sheet.  This  component  allows 
sporadic  merging  to  occur  where  the  merging  rate  is  controlled  by  the  magnetic 
noise  power  density  spectrum  and  the  value  of  B in  the  tail  lobes.  Upon  expansion 
of  the  plasma  sheet,  the  harmonic  component  cannot  be  maintained  and  merging 
stops.  This  process  is  repeated  periodically  until  enhanced  turbulence  damps  all 
coherent  (oscillatory)  motion. 

B-Field  Data:  The  two-component  steady-state  and  time -dependent  model  is 
supported  by  magnetic  field  data  taken  in  the  plasma  sheet  from  Explorer  34 


9.  Bowling,  S.  B.,  and  Wolf,  R.  A.  (1974)  The  motion  and  magnetic  structure  of 
the  plasma  sheet  near  30  R^,  Planet.  Space  Sci.  22:673. 
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(Bowling  and  Wolf  ).  Bowling  and  Wolf  showed  that  the  average  x-component  of 
the  tail  magnetic  field  varied  linearly  on  either  side  of  the  neutral  sheet.  In  con- 
trast, Garrett10  showed  the  significant  presence  of  ULF  magnetic  fluctuations  over 
short  time-scales  using  the  same  data.  In  this  paper  we  take  the  point  of  view 
that  the  averaged,  linear  B-field  acts  on  the  protons  mentioned  above.  This  leads 
to  the  harmonic  particle  component.  The  ULF  fluctuations  are  considered  to  be 
caused  either  by  the  sporadic  merging  as  described  in  Section  3 or  by  the  oscilla- 
tion modes  described  in  Section  4. 

Definit ions:  The  plasma  sheet  is  taken  to  be  the  region  in  the  earth's  geo- 
magnetic tail  that  lies  above  and  below  the  neutral  sheet  and  is  characterized  by 
a magnetic  depression  (Garrett10).  The  neutral  sheet  is  defined  to  be  the  z = 0 
plane  where  B reverses  sign.  For  purposes  of  this  paper,  we  take  the  plasma 
sheet  and  the  current-carrying  region  to  be  the  same. 

Coordinate  System  Used:  The  positive  x-axis  points  away  from  the  earth. 

The  positive  y-axis  points  toward  dawn  and  the  positive  z-axis  points  northward. 


2.  THE  PLASMA  SHEET 

The  plasma  sheet  is  considered  to  be  a dynamic,  turbulent  medium.  The 
dynamic  components  average  out  over  long  periods  of  time  to  steady-state  values. 
However,  over  shorter  time-scales  the  dynamic  behavior  has  profound  effects. 
First,  we  examine  the  steady-rstate  component  and  determine  the  time-scales  of 
interest.  Then  a macroscopic  time -dependent  component  is  examined  and  con- 
ditioned for  its  presence  deduced. 

The  Steady-State:  The  basic  steady-state  equations  are  (Spitzer11): 

j 

VP  = ^ X B (1) 

E + V/a  XB  = r,  j + ± VPj  (2) 

where  P is  the  particle  pressure,  r\  the  resistivity,  and  n the  particle  number 
density. 

The  current  density  required  to  maintain  (1)  can  be  obtained  simply  by  taking 
the  cross  product  of  B 


10.  Garrett,  H.  B.  (1973)  ULF  magnetic  fluctuations  in  the  plasma  sheet  as 

recorded  by  the  Explorer  34  satellite,  J.  Geophys.  Res.  78:3799. 

11.  Spitzer,  Lyman,  Jr.  (1956)  Physics  of  Fully  Ionized  Gases.  Interscience 

Publ.  — 


9 


— 
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= c 


B X VP 


(3) 


B 


We  assume  the  B-field  components  inside  the  plasma  sheet  to  be 


B = B = 0 

z y 


; 

and 


-2  B z 
2 


(4) 
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The  latter  is  based  on  the  measurements  of  Bowling  and  Wolf.  Here  Bq  is  the 


x magnetic  field  intensity  in  the  tail  lobes  and  I the  thickness  of  the  plasma  sheet 
in  the  z direction.  Equation  (1)  can  be  written  in  the  more  usual  form 


VP  = -V 


(*) 

Inserting  (5)  into  (3)  using  (4)  yields 


(5) 


cB 


2 n 


(6) 


which,  of  course,  can  be  directly  obtained  from  the  Biot-Savart  law. 
From  (l)  and  (2),  the  resistivity  is  near  the  neutral  sheet 


i ■ 

i 


n = Vjy 


m^ 


(7) 


ne 


Here  v is  the  collision  frequency  and  m.  instead  of  m is  used  because  protons 
carry  most  of  the  current  (Hones  et  al^). 


The  resistivity,  therefore,  is  a sensitive  measure  of  proton  scattering  by 

13  -3 

turbulence  (Akhiezer  et  al  ).  For  illustrations  we  choose  B = lOy,  n = 0.  3 cm  , 


1 = 2.3  B and  E ~ 10-5  volts/m.  Under  these  quiet-time  conditions 
e y ' 


12 


Hones,  Jr.,  E.W.,  Bame,  S.  J.,  and  Asbridge,  J.  R.  (1976)  Proton  flow 
measurements  in  the  magnetotail  plasma  sheet  made  with  Imp  6, 

J.  Geophys.  Res.  81:227. 

13.  Akhiezer,  A. I.,  and  Akhiezer,  I.  A.  (1975)  Plasma  Electrodynamics  Volume 
2:  Non-linear  Theory  and  Fluctuations.  Pergamon  Press,  Inc. , New  York. 


—I 


10 


r = \ /v  = 25  sec 

_3 

For  comparable  disturbed  conditions  we  take  Bq  ~ lOy,  n ~ 0.  3 cm  , 

I ~ 1500  km  and  E^.  ~ 2 X 10  ^ V/m.  In  this  case 

o 

j = 10.  6 na/m 

v ~ 0.  9 collisions/sec  (9) 

r = l/v  = 1.  1 6 3C 

Equations  (8)  and  (9)  give  time-scales  within  which  a coherent  time-dependent 
component  can  exist  before  being  dissipated  by  turbulence.  The  effective  electron- 

ion  collision  frequency  associated  with  ion-acoustic  turbulence  is  ~0.  4 T /T. 

14  e i 

(Biskamp  ) for  the  disturbed  conditions  given  above.  This  is  on  the  same  order 

as  (9). 

Time -Dependent  Behavior:  In  Appendix  A we  have  solved  the  Lorentz  equation 
exactly  in  the  limit  of  zero  electric  field.  This  limit  is  considered  appropriate 
since  E-field  acceleration  should  be  small  between  individual  collisions.  The 
exact  solutions  take  the  form  of  Jacobi  elliptic  functions  which  are  well  defined  in 
standard  texts  of  Mathematical  Physics.  As  shown  in  Appendix  A and  Figure  Al, 
the  topology  of  individual  particle  trajectories  is  dependent  on  the  numerical  value 
of  the  modulus  associated  with  the  Jacobi  functions.  The  particle  energy  (Ecp)  that 
corresponds  to  the  modulus -squared  equal  to  one  divides  particle  trajectories  into 
two  main  groups;  those  of  lower  energy  that  do  not  cross  the  neutral  sheet  and 
those  of  higher  energy  that  do.  The  critical  energy,  Ecp,  for  particles  that  have 
their  turning  points  on  the  plasma  sheet  boundaries  is  given  by 


cp 


m 2 .2 
128  wci 


(10) 


14.  Biskamp,  D.  (1973)  Collisionless  shock  waves  in  plasmas,  Nucl.  Fusion 
13:719.  " 
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where  - e BQ/mc  is  the  gyrofrequency  and  m the  particles  mass.  In  Appen- 
dix A it  is  shown  that  for  protons  Ecp  » kT  when  I ~ 2.3  Rg,  kT  = 2 keV  and  that 

E « k'T  when  i ~ 1500  km.  In  other  words,  the  Lorentz  equation  predicts  that 
cp 

the  majority  of  protons  will  follow  trajectories  as  shown  in  Figure  Aid.  rather 
that  those  in  Figure  Ala.  as  the  plasma  sheet  is  compressed.  The  trajectory 
shown  in  Figure  Aid.  is  harmonic  and  it  is  the  turning  on  and  off  of  this  harmonic 
motion  as  the  plasma  sheet  contracts  and  expands  that  regulates  substorm  activity. 

The  periodicity  of  the  north-south  proton  motion  is  directly  found  by  setting 
the  first  argument  in  the  Jacobi  function  cn  (see  Eq.  (A7))  to  4K.  (This  particular 
Jacobi  function  has,  by  definition,  a period  of  4K  in  the  first  argument.)  Here  K 
is  the  complete  elliptic  integral  of  the  first  kind.  The  result  is 


m 


2 2 

where  V is  the  proton's  velocity,  k = zm  uc/2V  i and  zm  = the  distance  of  the 

particle's  turning  point  from  the  neutral  sheet.  The  complete  elliptic  integral  of 

o 

the  first  kind  becomes  singular  as  k -»  1,  which  corresponds  to  the  case  shown 
in  Figure  A lb. 

Now  let  us  develop  an  equivalent  fluid  representation  of  the  single -particle 
motion  shown  in  Figure  Aid.  In  the  cold  plasma  limit  and  in  the  absence  of  bulk 
flow  the  fluid  equation  of  motion  is 


p ff  = kxS  = 


This  equation  implicitly  assumes  that  the  cold  protons  are  the  sole  source  of  B. 
For  the  field  specified  in  Eq.  (4),  Eq.  (12)  also  reduces  to  that  of  a harmonic 
oscillator 


p z = - 


where 


“ " f VA 


VA  = 


12 


The  solution  for  a fluid  element  starting  at  z = zm  with  zero  initial  velocity  in  the 
Z-direction  is 


z = z„  cos  ut 
m 


(15) 


If  Eq.  (14)  is  to  represent  the  same  dynamics  as  Eq.  (11),  then  a correction  term 
must  be  applied  to  account  for  the  singular  behavior  of  K.  This  correction  term 
is 


fry 

2K(k) 


(16) 


such  that  id'  -»u,  K -»»/2,  as  k -»  0.  The  case  k — 0 is  the  large  velocity  limit. 
The  condition  that  Eq.  (16)  is  an  alternate  fluid  representation  of  the  s ingle  - 
proton  harmonic  motion  shown  in  Eq.  (11)  is 


T'  = 2»/u'  = r 


which  reduces  to 


j = nev  = cBq/2  ifl 


(17) 


(18) 


Equation  (18)  may  be  stated  in  terms  of  a theorem.  If  the  current  density  in 
the  plasma  sheet  is  taken  as  nev  where  n is  the  proton  number  density  and  v the 
total  mean  velocity  in  the  Z - Y,  plane  then  the  time -dependent  cold  plasma  equa- 
tion of  motion  will  be  that  of  an  oscillator  with  a period  in  the  Z (north-south) 
direction  which  is  equivalent  to  that  obtained  by  an  exact  solution  of  the  Lorentz 
equation.  The  physical  picture  of  an  ensemble  of  protons  moving  with  random 

phase  throughout  the  plasma  sheet  can  be  replaced  by  a continuum  or  field  of  inde- 

1 5 

pendent  harmonic  oscillators  (Henley  and  Thirring  ) as  defined  by  the  cold  plasma 
Eq.  (13).  This  is  the  equivalent,  time -dependent  fluid  representation  of  the  single- 
particle motion  in  phase  space.  When  averaged  over  a surface  element,  both 
treatments  lead  to*the  same  result  since  they  have  identical  dynamical  properties. 

The  physical  significance  of  the  cold  plasma  approach  can  be  considered  as 
follows:  Introduce  a finite  collision  frequency  into  the  above  treatment.  Over  long 
time  periods  the  single  particle  motion  will  be  averaged  or  washed  out  by  collisions. 
Howevbr,  let  us  consider  shorter  time  scales.  The  number  of  expected  collisions 
over  a given  time  scale  follows  Poisson  statistics  so  that  the  number  of  protons 


that  do  not  collide  increases  exponentially  as  the  time  interval  decreases.  Even- 
tually, one  reaches  a limit  where  the  cold  plasma  equations  are  valid.  In  the 
present  case,  the  cold  equations  need  only  be  valid  over  the  time  interval  required 
for  a proton  to  cross  the  plasma  sheet  in  the  z-direction.  For  Bq  = 30 y, 

1 = 1500  km,  E = 5.0  keV  this  is  ~ 2 sec  (Eq.  (ID)  which  is  comparable  to  the 
estimated  mean  collision  time  of  ~ 1 sec  (Eq.  (9)). 

The  protons  in  this  picture  can  be  considered  test  particles  (Akhiezer  and 
Akhiezer13)  that  are  scattered  by  plasma  fluctuations  generated  by  collective 
phenomenon.  The  combined  effects  of  all  such  generated  phenomena  is  incorporated 
in  the  numerical  value  of  the  collision  frequency,  v.  A finite  v introduces  a 
Langevin  collision  term  on  the  right-hand  side  of  the  Boltzmann  equation  which 
implies  a net  momentum  transfer  from  protons  to  electrons  via  plasma  fluctua- 
tions. With  such  a term,  Eq.  (13)  becomes  the  well-known  equation  for  a damped 
oscillator.  Thus,  the  importance  of  single  particle  motion  in  controlling  the 
dynamics  of  the  tail  is  modulated  by  the  degree  of  plasma  turbulence  generated  by 

instabilities.  In  Section  4 we  show  that  for  critical  damping  (u  = v/2),  there  will 

2 

be  a rapid  heating  of  electrons  near  the  neutral  sheet  at  a rate  ~ (Bq/8»t)V^. 

We  now  derive  a simple  criteria  for  the  presence  of  the  harmonic  motion 
across  the  neutral  sheet.  Consider  the  particles  as  having  two  independent  oscil- 
latory modes:  one  being  the  usual  gyromotion  represented  by  u>c,  and  the  other 
being  the  north-south  harmonic  motion  represented  by  u>.  The  north-south  compo- 
nent will  dominate  when 

u)sHi=l/2ii)  (19) 

' c 

which  reduces  to 

nf2<8.  15  X 101 5 proton/cm  (20) 

3 

For  a nominal  density  of  0.  3 cm 
l s 1634  km 

As  in  the  single-particle  treatment,  a decrease  in  f is  required  to  trigger  the 
north-south  harmonic  motion.  Based  on  the  above,  we  use  ~1500  km  as  a nominal 
value  for  the  plasma  sheet  thickness  during  onsets. 

Equilibrium  Distribution  of  the  Harmonic  Component:  The  tir  e-dependent 
component,  because  of  Its  oscillator  properties,  can  lead  to  a steady-state  dis- 
tribution. For  any  oscillator  one  may  eliminate  the  independent  variable  t from 
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I 


the  (z,  z)  equations.  Then  an  observer  at  coordinate  point  z will  measure  a 
density  n (z,  z ) from  protons  originating  at  zm»  That  is: 


n (z. 


z ) = 
m 


const. 


const. 


u (z; 


2. 

z ) 


172 


(21) 


Now  only  protons  with  turning  points  greater  than  z can  pass  through  z.  There- 
fore, 


n'(z) 


i/2 

/ “ (z«  zm)  dzm 


n'  (z)  dz  = a ni  /2 


(22) 


where  u is  defined  in  (14).  This  leads  to  an  enhanced  density  at  z = 0 relative  to 
that  at  z = ±f  /2.  The  integration  of  n'  over  z from  0 to  i /2  determines  the  con- 

_3 

stant.  This  constant  is  equal  to  4onu/i^  where  n « 0.  3 cm  and  a is  defined  by 
(24). 

Multiple  Onsets:  Consider  a localized  impulse  imparted  to  the  external 
boundary  of  the  north  tail  lobe;  This  impulse  will  propagate  through  the  tail  lobe 
via  a compressional  wave  as  described  in  Section  4.  Upon  reaching  the  plasma 
sheet  boundary,  it  will  impart  additional  energy  to  protons  just  turning  at  the 
interface  (see  Figure  A Id. ).  These  protons  will,  if  not  scattered,  transfer  the 
impulse  to  the  opposite  tail  lobe  boundary.  The  presence  of  the  north-south  har- 
monic motion,  therefore,  will  dynamically  couple  the  two  lobes  and  the  entire  tail 
will  respond  as  a unit.  Any  slowing  down  or  scattering  of  the  coherent  protons 
would,  of  course,  prevent  them  reaching  the  opposite  boundary  and,  hence,  the 
coupling  would  be  broken. 

The  probability  that  a proton  will  not  suffer  a collision  during  a time  interval, 
t,  is  given  by  the  Poisson  distribution, 

p (0,  yt)  = e_I/t  (23) 

Returning  to  the  fluid  picture  where  the  harmonic  motion  is  represented  by  a 
field  or  continuum  of  independent  oscillators  the  fractional  number  of  these  oscil- 
lators that  are  undamped  during  a half-period  is  given  by  . The  effective 
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oscillator  strength  of  the  plasma  sheet  in  terms  of  ability  to  coherently  transfer 
information  from  one  boundary  to  the  other  is  therefore  given  by 

Kp  = u'2  e~*vlu'  =au'2  <24) 

The  parameter,  Kp,  is  the  equivalent  spring  constant  of  the  plasma  sheet.  In 
Section  4 an  equivalent  constant  will  be  derived  for  the  tail  lobes.  The  eigenmodes 
of  such  a system  are  then  examined  in  terms  of  multiple  onsets.  The  periodicity 
of  the  onsets  is  modulated  by  v through  Eq.  (24). 

But  first  in  Section  3 we  derive  a theory  for  merging  based  on  magnetic  insta- 
bilities that  are  generated  by  localized  currents.  These  localized  currents  are 
presumed  generated  by  particle  motion  such  as  shown  in  Figure  Alb.  This  motion 
is  also  significant  when  the  plasma  sheet  is  sufficiently  compressed. 


3.  MERGING 


In  the  previous  section,  it  was  shown  that  a transverse  harmonic  component 
can  exist  if  the  plasma  sheet  is  sufficiently  compressed.  In  the  present  section, 
we  derive  a theory  showing  it  is  possible  to  have  direct  field  line  merging  if  the 
magnetic  instabilities  are  large  enough.  This  direct  merging  is  advanced  as  an 
explanation  for  the  sporadic  reconnection  at  maximum  compression  reported  by 
Pytte  et  al.  1 Energy  conservation  is  also  discussed. 

What  is  meant  by  direct  merging  is  shown  in  Figure  1.  A localized  current 
fluctuation  causes  deformation  of  the  field  lines  on  the  plasma  sheet  boundaries. 

As  these  field  lines  convect  into  the  plasma  sheet  they  become  statically  unstable. 
It  is  argued  that  if  (19)  is  satisfied  then  the  motion  of  the  attached  electrons  can 
cause  these  field  lines  to  uniformly  and  smoothly  connect  across  the  x-axis  and 
merge  into  two  separate  field  lines.  Criteria  for  this  process  to  occur  is  derived. 

This  problem  is  approached  by  determining  consistency  conditions  on  the  de- 
formation itself.  The  corresponding  B-field  inside  the  plasma  sheet  is  then  de- 
rived. We  assume  that  in  the  vicinity  of  the  deformation  that  B is  a weak  function 
of  Z.  That  is: 


!fx  „ ffz 
9z  ~ dz 


~ 0 


which  implies  from  the  constraint  V • B = 0 that 
dB 

— L = 0 or  Bx  = const.  = -Bq 


(25) 


(26) 
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The  z-component  is  then  determined  directly  from  the  slope  of  the  deformation 


f 


! 

* 

\ 


■ 


t 


i 


B = -B  f'(x)  (27) 

z o 

Equations  (26)  and  (27)  represent  the  simplest  form  for  a B-field  line  next  to  that 
for  a uniform  field.  We  now  build  a theory  assuming  that  (26)  and  (27)  properly 
describe  the  deformations.  Equation  (27)  suggests  that  current  fluctuations  add 
B-field  energy  to  the  deformation  through  the  creation  of  a Bz  component. 


DIRECT  MERGING  BY  CURRENT  FLUCTUATIONS 


Figure  1.  A Schematic  Representation  of  the  Earth's  Magnetotail.  The  solid  lines 
show  the  tail  lobe  field  lines  being  deformed  by  a localized  interruption  of  the 
cross-tail  current  near  z = 0 (see  text).  The  dotted  lines  represent  the  connected 
field  lines  inside  the  plasma  sheet.  Merged  field  lines  are  represented  by  the 
dash-dot -dash  lines 


B-field  Lines  Inside  the  Plasma  Sheet.  The  stresses  acting  on  the  deforma- 
tion have  opposing  x-components  (Figure  1).  This  fact,  together  with 
Bx  (x  - 0)  < Bq  in  the  plasma  sh«:et,  implies  that  the  convected  field  lines  will  con- 
tinue to  deform  and  stretch.  We  assume  that  these  field  lines  will  be  closely 
approximated  by  a family  of  field  lines  that  (1)  join  smoothly  to  the  deformation 


(2)  are  divergence  free  and  (3)  are  symmetric  about  the  x-axis. 


(28) 


Bx  = -27^y  • '<*>=i0  + 2«x) 


where  i(x)  is  the  separation  between  the  plasma  sheet  boundaries  in  the  region  of 
the  deformation.  The  parameter  l denotes  the  minimum  separation  at  x = 0. 

The  z-component  of  B is  assumed  to  have  the  functional  form. 

Bz(x,  z)  = -BQf'(x)  g(x,  z)  (29) 

where  g(x,  z)  is  to  be  determined  subject  to  the  constraint  g(z  = t (x)/2)  = 1.  Using 
(28)  and  (29)  V-B  - 0 implies 

g,  z =-^|  (30) 

l 

which  upon  integrating  yields 


g(x,  z)  = + 1/2 

r 


The  factor  of  1/2  comes  from  the  boundary  condition  at  the  deformation.  The  B2 
component  is,  therefore. 


Bz(x,  z)  = - [2z2/i  2 + 1/2)  BQf'(x)  (32) 

The  point  is  that  as  z — 0,  Bx  — 0,  however,  Bz  does  not  vanish  at  z = 0 but 
has  a spectrum  of  values  depending  on  the  steepness  of  the  deformation.  It  is 
argued  that  the  convected  field  lines  originating  with  deformations  on  opposite 
sides  of  the  neutral  sheet  may  join  by  superposition  to  form  a new  family  of  field- 
lines; one  set  earthward  and  the  other  tailward  of  the  merging  point  (see  Figure  1). 

1 6 

Merged  Field  Lines;  The  simplest  form  for  a merged  field  (Eastwood  ) is 


Bx  = -2  Bq  z/f 


Bz  = i Bzl  = const. 

By 


16.  Eastwood,  J.W.  (1973)  Consistency  of  fields  and  particle  motion  in  the 
"Speiser"  model  of  the  current  sheet.  Planet.  Space  Sci.  20:1555. 


The  shape  of  the  field  lines  Is  determined  by 


That  is,  a parabolic  shape  about  the  x-axis.  Now  in  the  region  of  the  deformation 
an  analysis  similar  to  (31)  is  carried  out  except  now  the  boundary  condition  on  B 


where  Bz-  is  the  residual  value  of  Bz  (Bowling  and  Wolf,  and  Behannon 
leads  to 


Outside  the  deformed  region  where  f'(x)  = o Eq.  (36)  reduces  to  (33) 
At  some  point  x,  Eq.  (32)  will  be  indistinguishable  from  (36)  if 


In  this  case  the  B-field  components  of  a field  line  originating  on  the  deformations 

will  be  the  same  as  those  of  a merged  field  line.  Since  the  values  of  P varies 

from  o at  the  bottom  of  the  deformation  to  some  maximum  value,  P , at  the 

m 

inflection  point  then  (37)  can  be  expressed  as  the  inequality. 


Let  us  refer  to  Figure  2 for  the  physical  implications  of  this  result.  Assume 
there  exists  a power  density  spectrum  for  magnetic  fluctuations  as  a function  of 
P . This  is  represented  by  the  curve  on  the  left-hand  side  of  the  figure.  Quiet 


17.  Behannon,  K.  W.  (1970)  Geometry  of  the  geomagnetic  tail,  J.  Geophys 
74:743.  


(c) 


Figure  2.  Assumption:  There  Exists  a Power  Density  Spectrum  for  Magnetic 
Fluctuations  as  a Function  ffo.  (a)  A typical  magnetic  instability  (deformation). 
The  maximum  slope,  f^,  is  at  the  inflection  point;  (b)  A hypothetical  power  spec- 
trum as  a function  of  ^m;  (c)  An  idealized  model  showing  how  merging  on  the  day 
side  causes  B0  to  increase,  thereby  decreasing  the  threshold  at  which  merging 
occurs 


time  conditions  are  represented  by  point  B.  As  merging  increases  on  the  dayside, 

B„  increases  in  the  tail  lobes,  therefore,  decreasing  the  threshold  value  of  f'  at 
o m 

which  merging  can  take  place.  This  is  denoted  by  point  A.  In  this  manner  the 
nightside  merging  rate  is  controlled  by  the  dayside  rate.  The  process  is  very 
analogous  to  a variable  threshold  discriminator  where  the  threshold  is  regulated 
by  the  dayside  merging  rate  through  Bq.  An  increase  in  Bzi  would  decrease  the 
merging  rate,  thereby  allowing  the  temporary  build-up  of  magnetic  energy  in  the 
tail  lobes. 
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. , . cBo  L 4z2  fl2~|  cBo  T + 4z2~| 

]y<x. z)  = -w  |_1  + 7f  J |_  + 72- J 


At  x = 0,  z = 0,  f1  = 0 and  f"  (o)  is  equal  to  the  curvature  K of  the  field  line  at 
the  boundary  so  that 


jjo,  o)  = - 


The  current  density  is  positive  if  the  curvature  of  the  deformation  is  large  enough. 

K > 4/1  (41) 

Alternatively  the  presence  of  a positive  flowing  current  implies  positive  curvature. 
A positive  flowing  current  fluctuation  is  supplied  by  the  single  particle  motion 
described  in  Appendix  A and  in  Figure  Alb. 

Energy  Conservation.  The  acceleration  of  plasma  by  the  relaxation  of  the 
merged  field  lines  will  now  be  derived.  It  is  shown  that  the  magnetic  energy  is  con- 
verted into  plasma  kinetic  energy.  For  simplicity,  we  assume  Bz  = Bzi  through- 
out the  plasma  sheet  where  Bz  points  northward  earthward  of  the  merging  region 
and  southward  tailward  of  the  region.  The  Bx  component  is  taken  as  shown  in 
(33).  In  that  case,  it  was  shown  that  the  merged  field  lines  have  a parabolic  shape. 
The  orthogonal  plasma  trajectories  are  determined  by 


which  can  be  integrated  and  solved  for  z. 


si  exp  [-  b“T  (x  - xi>] 


The  accelerated  plasma  asymptotically  approaches  the  x-axis  as  x — 00. 


i 


Energy  conservation  requires. 


x,  z 


1/2  p (v2  - vf)  = I f <jyX  Bi- 


ds , 


(44) 


Xi>  zi 


where  i-subscripts  denote  initial  values  and  ds^  is  the  element  of  length  perpen- 
dicular to  B.  Now 


2 'I2 

B ■ (Bzi  * 4Bo  -fr) 


(45) 


(46) 


Inserting  the  defined  variables  from  (45)  and  (46)  into  (44),  one  finds 

,2 


l/2p(v2-vf)  = -|^ 


x,  z B2. 

— z/z- 
4?r  ' i 

xi*  zi 


(47) 


Using  (43),  we  find  the  second  term  on  the  R.  H.  S.  is  equal  to 

F JyBzi(x-xl)  (48) 

which  is  interpreted  as  work  performed  by  the  current  in  moving  the  particles 
from  (x.,  z.)  to  (x,  z).  Here  j is  the  unperturbed  current.  In  the  limit  x — «, 
z -*  o by  (43)  and  Bx  o (33).  Eq.  (47)  becomes 

1/2  p V2  = 1/2  p V2  + B2/8*  + I jy  Bzi(x  - x.)  (49) 

We  have  replaced  Bx  (x.,  z.)  by  Bq  since  we  have  shown  that  in  the  merging  region 
the  Bx  field  energy  becomes  converted  into  particle  kinetic  energy. 
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4.  MAGNETOTAIL  OSCILLATIONS 

Multiple  onsets  are  a distinctive  feature  of  magnetospheric  substorms  (Pytte 
et  al*'.  In  this  section,  a theory  is  derived  to  explain  the  observed  intervals 
between  onsets.  This  theory  predicts  a proton  collision  frequency  consistent  with 
(9).  Thermal  heating  of  the  plasma  sheet  is  also  discussed. 

It  was  shown  in  Section  2 that  a substorm  onset  could  lead  to  harmonic  par- 
ticle motion  in  the  plasma  sheet  near  the  region  of  maximum  compression.  We 
now  argue  that  this  motion  couples  to  the  tail  lobes  causing  the  tail  to  respond  as 
a complete  system. 

Compression  of  the  tail  magnetic  field  generates  transverse  compressional 
waves  across  the  tail  lobes  that  obey  the  equation  (Spitzer^1). 


v2e. 


' (vT *■?) 


92  E 


9t 


(50) 


92E. 


1 » 

~ ~2 2~ 

V2  9t2 


where  is  again  the  Alfven  velocity.  Now 


z = cEy/B  = =£-  E ei(kz-ut) 
J ' o B o 
o 


(51) 


so  that 


• ••  2 2. 

z = -lc  V . z 
A 

where  k is  the  wave  vector.  This  reduces  upon  integration  to 


• • 2 2 
z = -kz  VZ  z 


(52) 


(53) 


which  is  again  the  equation  for  a harmonic  oscillator.  The  reduction  of  the  wave 

1 5 

equation  to  that  of  a harmonic  oscillator  is  well  known  (Henley  and  Thirring  ). 
Equation  (53)  determines  the  equivalent  force  constant,  K^,  in  the  tail  lobes. 


ka  = k2  V! 


(54) 
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The  wave-vector,  k,  will  only  have  certain  discrete  values  if  the  tail  lobi 

18 

act  as  wave -guides  (McClay  and  Radoski  ).  If  rQ  is  the  height  of  the  tail  lob 
the  z -direction,  then  the  compress ional  wave  will  reflect  at  the  plasma  sheet 
boundary  and  return  to  the  outer  surface  after  a time,  t,  where 


where  N is  a positive  integer.  This  simplified  derivation  leads  to  results  identical 

18 

to  those  obtained  from  more  sophisticated  approaches  (McClay  and  Radoski  ). 
Inserting  (58)  into  (54),  we  have 


The  continual  compression  of  the  tail  lobes  by  the  solar  wind  ensures  that  the 
natural  frequencies  represented  (59)  will  always  be  present. 

Now  the  oscillator  strength  in  the  plasma  sheet  as  taken  from  Section  2 is 


18.  McClay,  J.F.,  and  Radoski,  H.R.  (1967)  Hydromagnetic  propagation  in  a 
theta-model  geomagnetic  tail,  J.  Geophys.  Res.  72:4525. 


The  parameter  a is  a direct  measure  of  the  turbulence  in  the  plasma  sheet  where 

a small  value  of  a implies  a highly  turbulent  medium  (that  is,  a high  collision 

frequency).  The  turbulent  component  does  not  transfer  coherent  momentum  from 

one  tail  lobe  to  another.  The  same  values  for  B , p and  r are  assumed  ; Nath 

o o 
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tail  lobes  so  that  one  has  three  harmonic  oscillators  strung  in  a line,  with  tne 
outer  two  oscillators  having  identical  force  constants.  Figure  3 shows  the  con- 


ceived arrangement.  This  system  is  now  solved  in  the  limit  of  small  amplitude 
oscillations  in  Appendix  B. 


Figure  3.  The  Magnetotail  Modeled  as  Three  Cou- 
pled Oscillators.  is  determined  from  the  eigen- 
modes  of  transverse  compress ional  waves  in  the 
lobes.  The  current-carrying  region  is  defined  by  t 


Figure  4 shows  the  resulting  eigenmodes.  The  first  two  modes  correspond  to 

translation  of  the  plasma  sheet  without  involving  any  compression  or  expansion. 

The  last  two  modes  do,  however,  involve  compression  and  expansion  of  the  plasma 
-3  -3 

sheet.  For  n =10  cm  , B = 10  y,  r = 20  R and  with  the  assumption  K„  « K. 

* o o e r A 

the  positive  root  of  (B8)  becomes  using  (59) 

Ug  ~ \ / 2 ~ 0. 24  N rad/sec  (61) 

or 

T3  = 26/N  sec  (62) 

19 

Garrett  noted  an  oscillatory-type  magnetic  fluctuation  in  the  plasma  sheet  with  a 
period  on  the  order  of  50  sec. 

19.  Garrett,  H.  B.  (1972)  ULF  Magnetic  Fluctuations  in  the  Plasma  Sheet  as 

Recorded  by  the  Explorer  34  Satellite,  Masters  Thesis,  Rice  University. 
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EIGENMOOES  OF  THE  COUPLED  OSCILLATOR  MODEL 


C.  d. 


Figure  4.  The  Four  Possible  Eigenmodes  of  the  Coupled  Oscillator  System  Shown 
in  Figure  3 


The  fourth  eigenmode  shown  in  Figure  4 is  the  most  interesting  since  it  cor- 
responds to  both  tail  lobes  moving  toward  z = 0.  It  is  this  mode  that  one  expects 
to  be  excited  by  the  initial  compression  of  the  plasma  sheet.  The  associated  eigen- 
frequency  is  the  negative  root  of  (B8)  which  goes  to  zero  in  the  limit  K,,  -»  0 and  is 

r 

a measure  of  the  plasma  sheet's  compressibility  compared  with  that  of  the  tail 
lobes. 


This  fact  can  be  seen  more  clearly,  if  one  solves  for  Kp,  in  terms  of  the  eigen- 
frequency 


(64) 


Now  the  average  period  between  onsets  as  given  by  Pytte  et 


al 


is  12.  5 min  or 


= 0. 008  rad/sec 

The  value  of  using  the  above  values  for  the  tail  lobe  parameters  is  0.  029  so 
that  the  conditions  is  well  satisfied.  This  simplifies  (63)  to 


KF  = Ul 


(65) 


Equating  this  result  to  (24)  allows  the  mean  proton  collision  frequency  to  be  solved 
in  terms  of  u ^ and  u'  (Eq.  (16)) 


P)_ 


(66) 


A Consistency  Test:  The  plasma  sheet  is  assumed  to  have  the  same  represen- 
tative values  as  used  for  Eq.  (9)  with  kT  = 2 kev.  This  makes  w'  = 0.49  rad/sec. 
Using  (64)  for  v takes  on  the  value  of 


v = 1.27  collisions/sec 


(67) 


from  (66)  which  is  in  excellent  agreement  with  (9).  Note  that  (67)  was  derived 
from  the  period  between  onsets  while  (9)  was  derived  from  the  resistivity  relation. 
This  indicates  an  overall  consistency  of  the  theory  as  described  in  Sections  2 and  4. 

Taking  v on  the  order  of  one  collision  per  second  then  a (see  Eq.  (16))  is  equal 
to  0.  1 percent.  Recall  that  a is  the  fraction  of  protons  that  execute  at  least  one 
period  before  being  diverted  by  collisions.  This  is  verification  of  the  original 
hypothesis  that  the  harmonic  component  is  a small  fraction  of  the  turbulent  plasma 
medium  but  plays  an  extremely  important  role  in  terms  of  the  plasma  sheet 
dynamics. 

Thermal  Heating  of  the  Plasma  Sheet:  We  have  seen  that  a very  small  fraction 
of  protons  contribute  to  the  coupling^of  the  plasma  sheet  to  the  tail  lobes.  What 
about  the  other  protons  that  are  more  rapidly  thermalized  ? 

It  is  well  known  that  for  a damped  oscillator  that  the  frequency  for  critical 
damping  is  given  by  = v/2.  Note  that  in  the  present  case  = 0.  53  rad/sec  and 
v ~ 1 collision/sec  which  quantitatively  consistent  with  critical  damping.  The  rate 
at  which  energy  is  being  thermalized  is,  therefore,  given  by 


dE 

dt 


- P 


(68) 
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per  unit  volume.  The  general  solution  for  the  critically  damped  case  is 


z = 0 at  t = 0 have  been  used.  Equation  (69) 


where  the  initial  conditions  z = z 
can  be  differentiated  to  obtain 


which  when  inserted  into  (68)  and  time-averaged  over  a quarter-cycle  (t  = »/2u  ) 
gives 


0 to  z - t / 2 is  included  by 

m ' v 


The  contribution  from  all  elements  from  z 
integration. 


Inserting  u = 2/1  V.  into  (72)  gives  the  final  result 


which  is  the  rate  of  thermal  heating  per  square  centimeter.  Multiplying  by  the 
length  and  width  of  the  tail  over  which  heating  is  taking  place  gives  the  total  power 
input  for  one  hemisphere.  The  damping  of  the  JX  B motion  by  collisions  heats  the 
plasma  sheet  at  a rate  which  is  comparable  to  that  expected  if  magnetic  energy  was 
being  annihilated  at  the  Alfven  speed. 

Consider  the  analogous  case  of  an  electric  motor  in  a viscuous  medium  where 
an  external  "battery”  supplies  the  power  to  drive  the  armature.  In  the  same 
manner  the  source  of  the  cross  tail  electric  field  supplies  the  power  that  heats  the 

plasma  sheet  by  the  above  process.  Note,  however,  that  the  rate  of  heating  goes 
3 

as  Bq  and,  therefore,  is  controlled  by  the  dayside  merging  rate.  Heating  by 
collisional  damping  by  the  J X B force  will  exceed  joule  heating  if  Bq  > 22  y and  the 
parameter  values  given  in  Eq.  (9)  are  retained. 


5.  SUMMARY  AND  CONCLUSIONS 


In  this  model,  we  found  that  when  the  plasma  sheet  is  compressed  harmonic 
motion  acr  oss  the  plasma  sheet  is  enhanced.  This  is  true  for  both  the  fluid  and 
particle  treatments.  This  harmonic  component  dynamically  couples  the  plasma 
sheet  to  the  tail  lobes  and  acts  as  a substorm  triggering  mechanism.  The  charac- 
teristic frequency  of  the  oscillation  mode  determines  the  period  between  onsets 
which  is  a sensitive  measure  of  the  mean  proton  collision  frequency,  v,  in  the 
plasma  sheet.  Consistency  was  found  between  the  collision  frequency  obtained 
from  the  onset  period  and  from  the  resistivity  equation. 

In  the  fluid  approximation  the  plasma  sheet  can  also  be  thermally  heated  under 
compression  in  a manner  analogous  to  that  of  a critically  damped  oscillator.  This 
occurs  when  I ~ 4 V . 

Current  fluctuations  that  arise  from  changes  in  single  particle  motion  are 
considered  to  drive  B-field  fluctuations.  These  B-field  fluctuations  can  merge  by 
superposition  if  the  corresponding  deformation  is  severe  enough  and  if  merging  is 
taking  place  on  the  dayside.  Subsequent  relaxation  of  the  merged  lines  accelerates 
particles  parallel  to  B. 

We  conclude  that  the  above  ideas  provide  a useful  and  viable  model  for  the 
onset  of  magnetospheric  substorms. 
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In  this  appendix,  single  particle  motion  in  the  plasma  sheet  is  examined, 
equations  of  motion  (Speiser*)  in  the  present  system  of  coordinates  is 


By  integration  and  substitution  (Al) 


These  hold  for  a linear  z-dependence  in  B, 
can  be  rewritten  (Speiser1)  as 


Speiser,  T.W.  (1965)  Particle  trajectories  in  model  current  sheets,  1,  Analytic 
solutions,  J.  Geophys.  Res.  70:4219. 


i I 


The  case  for  long  times  was  treated  extensively  by  Speiser.  In  this  limit 


2e2EB  t , 

• • o i 

z ~ - 5 z = -u  z 


im2c 


(A3) 


which  is  the  equation  for  an  oscillator  with  a time -dependent  frequency.  For  com- 
pleteness, the  numerical  value  of  w is 


u = 0.  50  v/t  rad/sec 


(A  4) 


where  BQ  = 10  y,  E = 2mv/m  and  l - 1500  km.  This  is  analogous  to  the  fluid 
result  in  Section  2 of  the  front  text. 

In  the  alternate  limit  of  zero  electric  field  the  limit  (A3)  does  not  exist  and 
(Al)  reduces  to 


2eB 


y = - 


f me 


z z 


(A  5) 


2eB 


l me 


z y 


By  defining 


k2  = 


Zm  eBo 


2 v/mc 


(A  6) 


(A 5)  can  be  directly  solved  in  terms  of  elliptic  integrals  and  functions  (Byrd  and 


Friedman  ) for  three  cases  of  interest. 


Case  I: 


k2<l 


(A  7) 


\ m / 

(t?  • k) 


2.  Byrd,  P.  F.  and  Friedman,  M.D.  (1954)  Handbook  of  Elli^Uc  Integrals  for 


Engineers  and  Physicists,  Springer-Verlag  Publ. , Berl 
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k2  > 1 


y = - zm  E 


K?) 


+ 2k  vt  - vt 


<A9) 


z = ± z dn 
m 


(***•*) 


The  expressions  cn  and  dn  are  Jacobi  elliptic  functions.  E(u)  is  the  incomplete 
elliptic  integral  of  the  second  kind.  The  constants  of  integration  were  chosen  so 
that  y and  t equal  zero  when  z = ± z 

The  z-coordinate  has  the  following  possible  values 


Case  I: 


z s z < z 
m m 


Case  II: 


z„  < z < 0 or  0 < z < z 
m m 


(A  10) 


Case  III: 


- z 


/7TT 
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Note  that  protons  can  cross  the  z = 0 axis  only  for  Case  I (k  < 1).  See  Figure  Al. 
The  parameter  k is  inversely  proportional  to  the  total  particle  velocity  v.  There- 
fore, Case  111  corresponds  to  low  energy  protons  that  gyrate  solely  on  one  side  of 
the  z = o axis  (Figure  Ala).  Case  II  corresponds  to  protons  that  start  at  Z = ±Zm 
and  asymptotically  approach  z = 0 as  t -»  » (Figure  Alb).  Note  that  in  this  case 
the  proton  is  traveling  in  the  +y  direction  at  the  total  velocity.  Case  I corresponds 
to  protons  starting  at  ± z with  sufficient  energy  to  cross  z = 0 (Figures  Ale  and 
Aid). 


POSSIBLE  PROTON  TRAJECTORIES  IN  THE  PLASMA  SHEET  FOR  PROTONS 
PASSING  THROUGH  Z„»750km  B0«IOr  ^-ISOOkm 


Figure  Al.  Various  Proton  Trajectories  that  Illustrate  the  Several  Regimes  of  Vy 
as  Shown  in  Figure  A2.  (a)  Simple  grad  B drift  at  low  energies;  (b)  Rapid  loss  at 
selected  proton  energies  that  cause  an  opposing  current  along  +y;  (c)  Zero  drift 
velocity;  (d)  Higher  energy  protons,. that  add  to  the  current  causing  a pinch  effect. 


The  critical  energy  above  which  protons  execute  harmonic  motion  about  z = C 
is  found  by  setting  (A6)  to  1.0.  For  protons  starting  at  z = ±i /2  this  condition 
is  simply 


eB  . 

v = ~m§  s IT  wc 


If  Bq  ~ 10y  then  u>c  ~ 1 rad/sec  and  the  critical  energy  is 

E = 1/2  mv2  = 8.  15  X 10'15  f 2 ev  (i  in  cm)  (A12) 

The  critical  energy  is  a sensitive  function  of  the  plasma  sheet  thickness,  l.  Take 
the  thermal  energy  in  the  plasma  sheet  to  be  a couple  of  keV  then  for  large 
I (i  ~ 2.3  Rg)  ECp  = 17.  5 keV  which  says  that  few,  if  any,  protons  cross  the 
neutral  sheet.  Now  compress  the  plasma  sheet  to  I ~ 1500  km.  The  critical 
energy  is  now  0.  18  keV  which  implies  a majority  of  protons  execute  harmonic 
motion  across  the  neutral  sheet.  This  is  clearly  a single-particle  counterpart  to 
the  substorm  triggering  mechanism  discussed  in  the  text. 

Further  information  can  be  obtained  by  defining  the  proton  drift  velocity  in  the 
y-direction  as 


v = i i m y/t 
y t - « 


For  the  three  cases  discussed 


Case  I: 


[ , 2E  (k)  A 

^ -toet  j 


(Al# 


Case  II: 


v = +v 

y 


Case  III: 


2k  - 1 - 


lIDl 
G)  J 


Here  E(k)  and  E(l/k)  are  complete  elliptic  integrals  of  the  second  kind.  K is  the 
complete  elliptic  integral  of  the  first  kind.  Note  that  at  k = 0.  909  the  drift  energy, 
E1,  is  zero.  The  drift  velocity  is  plotted  in  Figure  A2  for  protons  passing  close 
to  the  plasma  sheet -tail  lobe  interface,  (z  = t/2).  The  plasma  sheet  thickness, 
l,  has  been  set  equal  to  1500  km.  The  drift  velocity  profile  for  other  plasma  sheet 
thicknesses  can  be  found  directly  from  Figure  A2  by  recalling  that  the  critical 
energy  scales  as  t -squared.  See  Eq.  (A  12). 
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PROTON  DRIFT  VELOCITY  IN  THE  PLASMA  SHEET  AS  A FUNCTION 
OF  ENERGY 


Figure  A2.  Proton  Drift  Velocity  in  the  Plasma  Sheet  as  a Function  of  Energy. 
This  curve  scales  as  the  square  of  the  plasma  sheet  thickness.  Protons  that  pass 
through  zm  = 750  km  (see  Figure  Al)  have  been  specifically  considered 


Let  us  draw  a physical  scenario  of  a substorm  onset  based  on  the  above  re- 
marks. The  plasma  sheet  is  again  assumed  to  have  a temperature  of  a couple  keV. 
During  quiet  time  conditions  (I  - 2.3  Rg)  E'  is  ~25  keV  and  from  Figure  A2  and 
Figure  Ala  most  if  not  all,  the  protons  are  drifting  toward  dawn  (+y)  without 
crossing  the  neutral  sheet.  Now  as  the  plasma  sheet  is  compressed  E'  becomes 
comparable  to  the  thermal  proton  energy  and  some  of  these  protons  now  gyrate 
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1 

act •.  ss  z = 0 and  drift  towards  dusk  (-y)  (see  Figure  Aid).  The  protons  that  re- 
verse their  drift  direction  add  to  and,  therefore,  by  the  Biot-Savat  law  cause  f 
to  continue  decreasing.  A decrease  in  t lowers  E'  even  further  leading  to  an 
equilibrium  condition  in  which  a further  lowering  of  E'  causes  a negligible  change 
in  jy.  The  conclusion  is  that  the  plasma  sheet  is  unstable  to  this  pinch  effect  once 
E'  becomes  some  multiple  of  the  thermal  energy.  For  illustration  purposes  we 
set  E'  = 6 keV  and  solve  for  f in  Eq.  (A 6).  The  corresponding  thickness  is  found 
to  be  — 1.  1 R_. 

£j 

Figure  A 2 leads  to  one  other  final  comment.  The  positive  spike  in  the  drift 

velocity  corresponds  to  protons  that  having  started  at  z = ±i  /2  asymptotically 

approach  but  do  not  cross  the  positive  y axis  (Case  II  and  Figure  Alb).  Therefore, 

o 

particles  very  close  to  ECp  (k  = 1)  will  cause  a positive  current  fluctuation.  This 
effect  is  identified  as  the  source  of  the  magnetic  fluctuations  described  in  Section  3. 


Solution  of  the  Eigenmodes 


The  model  is  defined  in  Figure  3.  The  potential  energy  of  the  entire  system 
can  be  written  as  (Goldstein1) 


KA 

V = <Z2 


r )2  + — (z  - z - i)2  + —£■  (z  - z - r )2 
o'  2 ' 3 2 n 2 'z4  z3  ro' 
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Using  the  usual  definition 

ni  = zi  " zio  (B2) 

i = 1,4 

where  the  zero  subscript  denotes  equilibrium  values.  Expression  (Bl)  then 
simplifies  to 

V = 1/2  Ka  (n2  - rjj)2  + 1/2  K (r,3  - n2)2  + 1/2  (n4  - pg)2  . (B3) 


1.  Goldstein,  H.  (1957)  Classical  Mechanics,  Addison-Wesley  Publ.  Co.,  Inc., 
Beading,  MA. 


A43  = "A13 

Finally,  for  the  fourth  root  u2  = K + - (K2  + K^)1/2  one  has 

/ -K  + (K2  + K2)1/2  \ 

A£4  = y J A 14  = _A34 

A44  = ”A  14 

Each  of  the  eigenvectors  are  of  unit  magnitude  (Goldstein1),  that  is 


(B12) 


4 

I 

i=l 


This  implies, 


An  = X'2 
A12  = l>2 


K, 


‘13 


2 K2  + K2 


- k 


1 1/2 


K, 


‘14 


2 K2  +K2 


- K ^ K2  + K2 


Tl/2 
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The  four  eigenmodes  are  schematically  shown  in  Figure  4. 


(B13) 


(B14) 

(B15) 

(B16) 


(B17) 
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MATRIX  FOR  THE  EI6ENM0DES 


Ka  -w2  -Ka  0 

-Ka  Ka  + K-w2  -K 

0 -K 

0 0 -K, 


Ka+K-w2  -Ka 


Ka-cu2 


Figure  Bl.  Energy  Matrix  for  the  Configuration  Shown 
in  Figure  3. 


